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ABSTRACT 

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black 
hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution 
about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation 
with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting 
nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though 
observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on 
conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that 
another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very 
small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss 
of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian 
contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to 
be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing 
GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission 
will be even more challenging than believed. From the other side, this behaviour could in principle be used as 
a signature of mass segregation in galactic nuclei. 



1. MOTIVATION 

A stellar mass black hole or neutron star executes ~ 10 5-6 
orbits during the final year of inspiral toward a ~ 10 6 M Q su- 
permassive black hole (MBH). The large number of cycles 
implies that a phase-coherent measurement of the inspiral, 
achievable through detection of low frequency gravitational 
waves, would be a tremen dously powerful probe of the space- 
time near a black hole dAmaro-Seoane et all 120071; [Hughes 
2009). Among other things, it would enable a precise determi- 
nation of the spin of the supermassive black hole, and a test of 
General Relativity that is independent of current constraints 
derived from pulsar timing data. 

There is no foreseeable instrument sensitive enough to de- 
tect gravitational waves from extreme mass ratio inspirals 
(EMRIs) over time scales comparable to the orbital period. As 
a consequence, realizing the astrophysical and gravitational 
physics promise of EMRIs requires an assurance that the in- 
spiral can be accurately modeled over many orbits using tem- 
plates calculated by solving the 2-body prob lem in General 
Relativity (for a review, see e.g., Barack 2009). It is therefore 
necessary to assess whether gas, stars or other compact ob- 
jects in the vicinity, could significantly perturb EMRI trajec- 
tories. In the case of gas, perturbations to stellar mass black 
holes or neutron starfl are securely negligible provided that 
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accretion on to the bla ck hole occurs in a low density, radia- 
tively inefficient flow (Narayan 2000). Such flows are much 
more common than dense accretion dis cs, which would yi eld 
observable phase shifts during inspiral dKocsis et al.ll201 ll) . at 
least at the relatively low redshifts where EMRIs may be ob- 
served. 

In this Letter, we quantify the nature and strength of possi- 
ble perturbations from point mass perturbers: low mass stars 
or compact objects in tight orbits around the supermassive 
black hole. Any perturbers are unlikely to orbit close enough 
to the EMRI to undergo strong interactions, so the regime of 
interest is one where the third body is relatively distant and the 
interaction weak. The Newtonian analog of this problem has 
been studied extensively in the context both of Solar System 
satellite evoluti on, and for transit timing variations of extra- 
solar planets (iDermott. Malhotra & Murravlll988t lAgol et al.l 
120051; iHolman & Murray! 120051; IVeras. Ford & Pavndl201 ll) . 
In Newtonian gravity, perturbations are strong only at the lo- 
cation of mean motion resonances, and these have the effect 
of inducing small jumps in eccentricity upon divergent res- 
onance crossing. This would already be interesting for the 
EMRI problem, since the jumps in eccentricity would result 
in a perturbation to the gravitational wave decay rate, and an 
eventual dephasing of the waveform. However, as we will see, 
the inclusion of post-Newtonian corrections changes the evo- 
lution qualitatively. Computing trajectories that include the 
two first-order non-dissipative post-Newtonian corrections, 
we find evidence of dependence on initial conditions in the 
evolution of the perturbed inner binary, such that arbitrarily 
small variations in the initial orbit lead to significantly differ- 
ent future behaviour. 

2. ASTROPHYSICAL LIMITS ON PERTURBERS 

Is it likely that a star or compact object will be present close 
enough to perturb the orbit of an EMRI? Excluding low mass 
MBHs (M. < 10 6 M Q ), where the stellar tidal disruption limit 
comes into play, the existence of perturbers is not excluded by 
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FIG. 1. — Estimates for the semi-major axis of the innermost perturbing 
body around a massive black hole, scaled to the hole's gravitational radius 
R g = GM,/c 2 . The red lines show the location of the innermost star, esti- 
mated assuming that stars of mass 0.3 Mq follow a single power-law cusp of 
index 7 in a galaxy on the M, -a relation. The green line shows the tidal dis- 
ruption limit for such stars. The blue lines show the average (upper) and 1% 
probability (lower) location of the next nearest EMRI, assuming uncorrected 
inspirals at a rate of 10~ 6 yr~' . 

elementary arguments. Neither, however, is it easy to calcu- 
late the probability distribution of perturbers, whose proxim- 
ity will depend upon the details of discreteness and relativis- 
tic effects very close to the MBH, and mass segr egation and 
EMRI injection mechanisms in galactic nuclei (fPfeto et alJ 
l2004t lFreitag"eTaIl l2006: Am aro-Seoane et alJl2004l) ~ 

Rather than face these difficulties, we limit ourselves here 
to order of magnitude estimates for the likely location of the 
nearest star and compact object. For stars, assumed to be of a 
single mass M*, we assume a cusp-like distribution with den- 
sity profile p oc R~~<, extending from the MBH to its radius of 
influence Rbh = GM,/a 2 . Here a is the velocity dispersion of 
the galaxy. Using the fact that the enclosed mass, M(R) ~ M. 
at R = Rbh, we find that the expected radius of the innermost 
star, is, 

!/(3-7) , „ , 2 
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(1) 



where R g = GM,/c 2 . This formula yields an explicit esti- 
mate for R\ once w e adopt a relation between M. and a 
dGultekin et al.l [2009h . For the location of the next nearest 
compact object (or EMRI), we use an even simpler approach. 
We calculate the expected semi-major axis for unc orrected 
inspirals due to gravitational radiation (Peters! [1964) . assum- 
ing near-circula r orbits an d rate /Vemri- Finally, we plot the 
tidal limit (e.g. lReesl[l988h for 0.3 M main-sequence stars. 

Figure Q] shows these estimates as a function of M # . For a 
standard cusp slope 7 = 1.75, there is likely to be a low mass 
stellar perturber within a few hundred R g for M, > 10 6 Mq. 
Similarly, if the EMRI rate is as high as 10~ 6 yr" 1 , there 
is a significant chance (at least a few percent) that a sec- 
ond compact object will be present between 10- 10 2 R g for 



10 6 Mq < M. < 10 7 Mq. Clearly, these crude estimates 
do not demonstrate that most EMRIs will be perturbed by 
third bodies, but they do suggest that perturbers may be close 
enough in some galaxies to motivate detailed consideration of 
their dynamical effects. 

3. METHODS 

We are interested in the secular effect of a star acting on 
an EMRI which will describe thousands of orbits in the de- 
tector bandwidth and slowly decay. The kind of effects on 
the wave that we are looking at are tiny, though detectable, 
and the mass difference between the two binaries (the MBH- 
EMRI and the MBH-star systems) is huge. We need there- 
fore a numerical tool capable of integrating the plunging orbit 
of the EMRI while inducing a minimal error in the integra- 
tion, since data analysis techniques can detect e.g. eccentric- 
ity di f ferences of the order A e ~ 10~ 3 dAmaro-Seoane et alJ 
l2010t iPorter & Sesanal l20l(it iKev & Cornish! 1201 lh. We 



hence have chosen to use a direct A^-body approach ( Aarsethl 
119991 12003h . the planet code, written by AarsetrQ. This 
is the most expensive method because it involves integrat- 
ing all gravitational forces for all three bodies at every time 
step, without making any a priori assumptions about the sys- 
tem. Our approach employs the improved Hermite integration 
scheme, which requires computation of not only the acceler- 
ations but also their time derivatives. Since we are simply in- 
tegrating Newton's equations directly, all gravitational effects 
are included. For the purpose of our study, nonetheless, we 
have included relativistic corrections to the Newtonian forces 
(the forces can be found in the same page in the toy cod^H). 
This wa s first implement ed in a direct-summation YV-body 
code by Kupi et al. (2006). For this, one has to add pertur- 
bations in the integration, so that the forces are modified by 



Newt periapsis shift energy loss neglected 

F = / ^+c' 2 F 2 +c~ A F i + c" 5 F 5 +0(c^) 



(2) 



1PN 



2PN 



In the last equation "PN" stands for post-Newtonian. We note 
that the perturba tions do not nee d to be small compared to the 
two-body force (Mikkola 199 71). The expres sions for Fi, F\ 
and F5 can be found in Blanche t & Fayd (1200 lh . their equation 
7.16. 

4. DISSIPATION OF ENERGY AND RESONANCES 

We first analyse the system by contemplating only the rel- 
ativistic effect of dissipation of energy; i.e. our simulations 
only incorporate the 2.5 PN correction term. We stop the in- 
tegration when the separation between the stellar BH and the 
MBH is a, = 5Rs c hw, which approximately corresponds to the 
limit where the PN approximation is not valid anymore. The 
inspiral down to this distance takes typically in our simula- 
tions some 440,000 orbits. 

In Fig.© the test stellar black hole of mass m, = 10M Q has 
been initially set in such an orbit that it is totally embedded in 
a LISA-like detector band (i.e. with an orbital period < 10 5 
sees, namely P m = 6 x 10 3 sees) and is hence an EMRI; its 
initial semi-major axis is a..; ~ 1.45 x 10~ 6 pc and its eccen- 
tricity e. ; = 0.05. The perturber, a star of mass = 10M Q is 

7 who, as is his admirable custom, has made the code publicly available 

|http : / /www . ast . cam . ac . uk/~sver re /web/pages /nbody ■ htm| 



I ftp : //ftp . ast . cam, ac ■ uk/pub/sverre/toy/README| 
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FIG. 2. — Upper panel: Results for the fiducial case using the direct- 
summation N— body integrator. The mass of the MBH is M, = 10 6 Mq, 
the mass of the stellar black hole is m. = IOMq . See text for more details. 
Lower panel: Same configuration but with an initial inclination of the star 
of u = 45° instead of 30°, i.e. the inclination triggers the Kozai mechanism, 
since u > 39.2° and the orbit is prograde. As mentioned in the previous case, 
even if the changes in eccentricity cannot be directly seen in the curve, they 
are of the order Ac, ~ 10~ 3 . 

initially on an orbit in which the semi-major axis has the value 
a+ i ~ 4.1 x 10~ 6 pc and the eccentricity at T = is e* = 0.5. 
The inclination of the system EMRI - star was set to 30° ini- 
tially in the upper panel. This constitutes our reference sys- 
tem. 

In the figure, the straight lines mark the condition P±/P,= 
A, with A an integer, the period of the star around the MBH 
and P, the period of the EMRI around the MBH; i.e. where 
the resonances occur. The first three resonances have an im- 
pact on e, which can be seen on the plot; later resonances do 
also affect e., with Ae, ~ 10~ 3 . We also note that in the upper 
panel one can see in-between smaller jumps; they correspond 
to higher-order resonances, P+/P, = 5.5, 6.5 and 7.5. 

We made the choice for an initial inclination of 30° to avoid 
another effect that introduces a change in the eccentricity. In 
the lower panel we have exactly the same system but for z* = 
45°. With this value, and the fact that the o rbit is progra de, the 
Kozai oscillation of eccentricity is present (Kozai 1962). Even 
if the eccentricity of the EMRI e, suffers the characteristic 
Kozai oscillations, the loci for the resonances still fulfill the 
condition P+/P, = integer. 

5. DOES THE FLAP OF THE STAR AT APOAPSIS SET OFF A 
TORNADO AT PERIAPSIS? 

In this subsection we address numerically the effect of in- 
cluding the relativistic periapsis shift along with the dissipa- 
tion of energy; i.e. the set of corrections as specified in Eq.©. 
As we show below, the effect of the periapsis shift changes 
completely the evolution of the system. In Fig.® we show 
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FIG. 3. — Fiducial case with energy dissipation and periapsis shift correct- 
ing terms for different initial inclinations of the perturber. The solid (red) 
curve corresponds to u = 30°, the long-dashed (green) to u = 30.001°, the 
short-dashed (blue) corresponds to the fiducial case plus a billionth of a de- 
gree, u = 30.0000000001° and the dotted (magenta) to the reference plus a 
Ur 13 of a degree, U = 30.0000000000001°. 

four cases. One of them corresponds to the reference system 
but taking into account the periapsis shift. We only display 
these examples but note that the behaviour is also chaotiqj 
for other nearby choices of i+. When using an initial inclina- 
tion of /* = 45°, which corresponds to the same situation as 
in the lower panel of Fig.© but taking into account the pe- 
riapsis shift, along with another case which is identical but 
for z* = 45.0000000001°, we find also a chaotic result which 
moreover eliminates the secular Kozai oscillation of e. 

We have systematically studied this chaotic behaviour by 
running hundreds of simulations in which we methodically 
increase in minimal differences an initial dynamical orbital 
parameter such as the inclination, semi-major axis or eccen- 
tricity. In all cases and parameters the evolution corroborates 
the chaotic behaviour of the system. We have also tested a 
mass for the perturbing star of 5 and 1.44 Mq, as well as dif- 
ferent values of e* (0.1, 0.3, 0.7 and 0.9), with similar results. 

In order to fence in the region within which the system is 
chaotic, we systematically increase the semi-major axis of the 
star and run the same experiment. We start with the same dif- 
ference in inclination at a slightly larger semi-major axis, and 
then regularly increase it until we reach one order of mag- 
nitude over the fiducial case, as we depict in Fig.©. The 
chaotic behaviour ceases at about one order of magnitude of 
the initial value of a+ in the reference case. 

6. QUANTIFYING THE DEPENDENCE ON INITIAL CONDITIONS 
OF THE SYSTEM 

In this section we present a way of characterizing the rate of 
separation of infinitesimally close trajectories systematically. 
To achieve this we compare our fiducial model with another 
case in which we set up the EMRI in an (almost) impercepti- 
bly different initial orbit (the initial difference is 2 x 10~ 10 pc, 
while the objects are moving on the same ellipse) and keep the 
same initial conditions of the MBH and the perturber. Hence 
EMRI in the second case differs only from the reference case 
slightly and has an initial distance separation of ro. We say 
that the two models are in phase provided that 

9 When we use the word, we do not follow the rigorous mathematical 
definition of chaos. We mean a strong dependence on the initial conditions. 
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FIG. 4. — Same as in Fig. (3) but we set initially the perturber at a larger 
and larger initial semi-major axis. From the top to bottom and from the left 
to the right, the semi-major axis of the perturber is a* = 4 X 10~ 6 pc, 6 X 10~ 6 
pc, 9 X 1(T 6 pc and 4.07243 X 10" 5 pc. Solid lines correspond to i* = 30° 
and dashed lines to U = 30.0000000001°. 



be at such a radius from the EMRI system that it will signif- 
icantly perturb it. From the standpoint of detection and data 
analysis, this is yet another complication of the problem and 
could even lead to the misinterpretation that nature's GR is not 
what we believe it to be. On the other hand, from the point of 
view of stellar dynamics, the detection of one of these systems 
would shed light on our current understanding of galactic dy- 
namics in general and mass segregation in particular. 
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r « r . (3) 
If the two different realizations reach a separation 

r^2a. (4) 

the EMRI bodies are moving out of phase, on entirely un- 
related orbits. We thence are able to estimate a characteris- 
tic timescale Td ep h for the system to become out of phase. In 
Fig.© we display the separation of the two systems for differ- 
ent distances to the perturber. From these figures we can mea- 
sure the value of a characteristic timescale Td ep h for a given 
a*. 

From the data points obtained in the upper panels of Fig.© 
we can then derive the relation displayed in the lower panel. 
For large enough distances, of the order of ~ 10~ 5 pc the two 
timescales converge and the system becomes deterministic. 

7. CONCLUSIONS 

In this paper we have addressed the role of a perturbation 
on an EMRI by a nearby star. The system depends extremely 
on minimal changes in the initial conditions (as small as a 
10~ 9 part in the inclination) lead to a very different dynami- 
cal evolution. In all cases, however, the Kozai mechanism is 
washed out by the periapsis shift, as o ne can expect (see e.g. 
iHolmanetaLl 119971 iBlaes et al.ll2002l) . For distances of the 
order of a* ~ 10~ 5 pc the system enters the chaotic regime, 
for perturbing masses as small as 1.44M©. While we can- 
not state clearly whether this will be a common feature for 
EMRIs, since the different dynamical and relativistic phe- 
nomena involved in the problem are many and not straight- 
forward (see for a rev iew lAmaro-Seoane et al.ll2007l and also 
lAmaro-Seoane 1201 ll for a dedicated review of the dynamics), 
it seems plausible that for a Milky Way-like galaxy a star can 
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FIG. 5. — Upper panels: From the left to the right and from the top to the 
bottom we show the separation r for a increasing separation of the perturbing 
star of 3.5 X 10~ 6 , 3.9 X 10~ 6 , 4.375 X 10~ 6 , and 4.5 X 10" 6 pc. The dashed 
line shows the critical distance 2a,. Note the different timescales in the lower 
panels. Lower panel: Tjj.pi, against distance to the perturber normalized to the 
gravitational radiation timescale of the isolated system Ti nsp ; i.e. the merger 
timescale without the perturber acting onto the binary MBH-EMRI. 
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